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A new Super-Kamiokande (SK) search for Supernova Relic Neutrinos (SRNs) was conducted using 
2853 live days of data. Sensitivity is now greatly improved compared to the 2003 SK result, which 
placed a flux limit near many theoretical predictions. This more detailed analysis includes a variety 
of improvements such as increased efficiency, a lower energy threshold, and an expanded data set. 
New combined upper limits on SRN flux are between 2.8 and 3.0 v e cm -2 s _1 > 16 MeV total 
positron energy (17.3 MeV Ej,). 



I. INTRODUCTION AND MOTIVATION supernova event is enormous. In a galaxy such as the 

Milky Way, it is estimated a supernova will occur around 

A. The SRN signal 



With supernovae each releasing on the order of 10 J 
of energy, 99% as neutrinos, the neutrino flux from each 
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2 or 3 times a century Although we would not see a 
significant neutrino burst from a supernova farther away 
than our galaxy and its satellites, all of the core collapse 
supernovae that have exploded throughout the history of 
the universe have released neutrinos, which, in the ab- 
sence of unexpected physics, should still be in existence. 
These neutrinos are herein called the supernova relic neu- 
trino (SRN) signal (or simply relics, though care should 
be taken, since this terminology is also sometimes used 
for big bang relic neutrinos), also sometimes referred to 
as the Diffuse Supernova Neutrino Background (DSNB). 

Many astrophysicists have considered the SRN signal, 
constructing models that predict both the flux and the 
spectrum. The first models were crafted even before 
SN1987A @-[4|], then after SN1987A with increasing in- 
terest and sophistication, with models such as: Totani et 
al.'s constant SN rate model @; Malaney's cosmic gas 
infall model @h Woosley and Hartmann's chemical evo- 
lution model [7j; Kaplinghat et al.'s heavy metal abun- 
dance model Q; Ando et al.'s LMA model Q Lunar- 
dini's failed supernova model [ll[ [Hj]; and the variable 
neutrino temperature formulation of Horiuchi et al. (6 
MeV and sometimes 4 MeV cases considered) [l3[ . Some 
examples of theoretical spectra are shown in Fig. [T] Al- 
though the normalizations vary, the general shape and 
slope of the predicted spectra are relatively similar. 

The SRN signal has never been seen. A paper was 
published in 2003 detailing the first search for the SRN 
events at SK pjj]. The basic method was to eliminate 
as many backgrounds as possible, then attempt to model 
the spectrum of the remaining backgrounds. A x 2 fit 
was performed on the energy spectrum of the final data 
sample, with two background components and one signal 
component. From this fit, a final model independent flux 
limit of 1.2 D e enr 2 s _1 for v e energy > 19.3 MeV was 
extracted, which was ^100 times more stringent than 
the previous world's best limit (l5j . This study is an up- 
date to that 2003 result, with more livetime, lower energy 
threshold, and significantly increased sensitivity. Other 
experiments have also produced SRN limits, such as SNO 
[ia ] and KamLAND [nj]. Recent models predict a SRN 
flux that is on the cusp of discovery [l3|, [l8| , motivating 
further efforts. 



B. The SK detector 

The Super-Kamiokande detector is a 50 kton water 
Cherenkov detector located in the Kamioka mine in 
Japan. The inner detector (32.5 ktons, 22.5 ktons fidu- 
cial) has an inner surface lined with ~11,100 (^5,200 for 
SK-II) 50 cm Hamamatsu PMTs. More information on 
the hardware and experimental setup of the detector, in- 
cluding calibration details, is provided by the following 
references: [l^-[2l|. 

The SK detector began data taking in April 1996, and 
was shut down for maintainance in July 2001. This pe- 
riod is referred to as SK-I (1497 days livetime). While 
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FIG. 1. Examples of theoretical SRN spectra. Flux is v e only. 
For heavy metal abundance model, average metal yield of one 
solar mass is assumed. 



refilling after the maintainance, on November 12, 2001, 
a major accident [22[ destroyed 60% of the PMT tubes. 
The surviving tubes were redistributed to cover the en- 
tire detector while replacement tubes were manufactured, 
and data taking was resumed in December 2002, running 
until October 2005 with reduced cathode coverage. This 
is the SK-II period (794 days). After the replacement 
tubes were installed, and the detector was back up to 
its nominal 40% coverage, data taking again continued 
from July 2006 until August 2008, marking the SK-III 
period (562 days). Finally, the detector electronics were 
upgraded [H, [24| , and data taking since September 2008 
is referred to as SK-IV. This study considers SK-I, SK-II, 
and SK-III data. 

The SK detector has many different triggers. When 
enough PMTs observe light within 200 ns (approximately 
the longest time it can take light to traverse the inner 
detector) a detector trigger occurs. Inner detector (ID) 
electronics triggers include the High Energy (HE) trigger 
(about 33 PMTs), Low Energy (LE) trigger (about 29 
PMTs, 100% efficient at 6.5 MeV), and Super Low En- 
ergy (SLE) triggers (variable threshold between 17 and 24 
PMTs). SLE events have the lowest threshold (which has 
changed over time) , and occur orders of magnitude more 
often than LE triggered events (~ 1 KHz as opposed to 
~ 10 Hz), requiring dedicated hardware for filtering and 
reducing the SLE data flow to manageable levels. The 
SLE trigger enables SK to see the lowest energy events. 
The outer detector (OD) also triggers if 19 or more OD 
PMTs fire within 200 ns. 

SK primarily sees SRN events via inverse beta decay 
{p e + p —> n + e + ). The next most visible mode for 
SRN interactions is about two orders of magnitude less 
frequent. The inverse beta decay cross section used is 
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described in references (2|| H(| . 



II. DATA REDUCTION 
A. Cuts 

Each SK PMT fires with a dark rate of -3-5 kHz; cos- 
mic ray muons penetrate the detector at ~2 Hz, and ^25 
atmospheric and solar neutrino events are identified ev- 
ery day. Many of these act as backgrounds to the SRN 
signal. With an expected rate of only a couple of SRN 
events a year, careful background reduction is necessary. 
Most backgrounds are removable, but a few remain af- 
ter all cuts and must be understood. The events that 
can be removed, and the cuts used to remove them, are 
described in this section. 

The software tools used for reconstructing event ver- 
tex, direction, and energy information are the same as 
used in the long-standing SK solar neutrino analysis 

(MM)- 



Noise reduction 

Many noise events can be separated from the main 
sample easily. These preliminary cuts are the same as 
those used for the long-studied SK solar analysis [13] > 
and are collectively called 'noise reduction'. 

Calibration events are labeled as such and are imme- 
diately removed. Events which trigger the OD are also 
removed, as these events have a charged particle enter- 
ing the detector from outside, while relic events should 
be fully contained. Electronics noise events are also cut 
using timing and charge considerations, with a signal ef- 
ficiency of almost 100%. 

The total charge deposited in the detector by the event 
is a quick criterion useful for separating out muons. 
Muons are in general much brighter than inverse beta 
decay positrons, and deposit much more charge in the 
detector. Any event with more than 800 photo electrons 
is removed (400 p.c. for SK-II), quickly eliminating most 
of the muons. 



Fiducial volume cut 

The SK ID is filled with 32.5 ktons of pure water. Ra- 
dioactive events occur in the detector from sources such 
as the surrounding rock wall and radioactive isotopes in 
the PMT glass itself. Most of these events are eliminated 
by definition of a fiducial volume (FV). In the relic anal- 
ysis, we have implemented the standard SK FV cut of 
removing the volume within 2 meters of the walls, leav- 
ing a FV mass of 22.5 ktons. 



Spallation cut 

The analysis lower energy threshold is determined by 
spallation, which is the process of a nucleus emitting nu- 
cleons when struck by a highly energetic particle. Cosmic 
ray muons interacting with an oxygen nucleus cause ra- 
dioactive spallation which mimics a real SRN signal. At 
lower energies, the unstable nuclei tend to have longer 
lifetimes, making them harder to tag through muon cor- 
relation. Eventually tagging becomes too inefficient for 
the remaining data to be useful, thus determining the 
lower energy threshold of the analysis. Spallation con- 
stitutes the bulk of the SRN background after the noise 
reduction cut. Elimination of spallation events is com- 
plex and quite costly in terms of signal efficiency, making 
it perhaps the most important removable background. 

To remove spallation, we correlate relic candidates to 
preceding muons. If the muon is responsible for the event, 
then the event should be located near the muon track, 
both in location and time. The main spallation cut uti- 
lizes a likelihood method: 

First, the muon tracks are reconstructed and catego- 
rized as: (i) single muons that travel through the detec- 
tor (~ 83% of muon events); (ii) multiple muons bun- 
dles (~ 8%); (iii) muons that enter from outside but 
stop within the ID (stopping muons, ~ 5%); and (iv) 
short track length muons, usually just clipping the cor- 
ner of the detector (~ 4%). Each of these muon types 
was found to produce spallation differently, and are thus 
treated independently. 

A dE / dx profile of each muon track is constructed by 
projecting back the amount of light seen by each PMT 
and determining where along the muon track it originated 
according to the measured arrival time of the light at 
the PMT. This profile of charge emitted along the muon 
track tends to peak for spallation producing muons ap- 
proximately at the position along the track the spallation 
occurred. Although the exact mechanism responsible for 
this peak is not known, the correlation is unmistakable 
in our data. 

Four spallation tagging variables are considered. The 
first spallation variable, At, is simply the amount of time 
by which the muon preceded the relic candidate. The 
second variable, Ltrans, is the transverse, or perpen- 
dicular, distance from the muon track to the relic can- 
didate (see Fig. [5]). The center of the nine consecutive 
bins (4.5 m distance) with the largest combined charge 
in the dE/dx histogram corresponds to a position along 
the muon track where the spallation is expected to have 
occurred. The distance along the muon track from where 
the spallation is expected to occur to where the relic can- 
didate exists is the third likelihood variable, longitudi- 
nal distance £long- Lastly, the value of the combined 
charge in the largest nine consecutive bins of the dE / dx 
histogram is our fourth likelihood variable, Qpeak- The 
more charge in the peak, the more likely the muon is to 
be making spallation. 

Probability density functions (PDFs) are formed for 
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FIG. 2. 

ables. 



Schematic explanation of spallation distance vari- 



each of the four spallation variables. For each muon 
categorization, the relic candidates were correlated to 
muons preceding the candidate in time (the 'data sam- 
ple'), and correlations were examined. A random sample 
was formed by taking the same correlations to muons 
immediately following the relic candidate in time. The 
random sample histograms were subtracted from the 
data histograms (yielding a 'spallation sample'), for each 
muon categorization. These profiles were parameterized, 
resulting in functions representing spallation (from the 
spallation sample) and accidental correlation background 
(from the random sample). These parameterizations, 
once normalized, are the PDFs, which are multiplied to- 
gether to give the likelihood (for example, see Fig. [3]). 

For each muon categorization, a cut on the likelihood 
was instituted. Cut values were tuned until no statisti- 
cally significant difference existed in the distributions of 
the data remaining after the cut compared to the random 
sample for the At and £trans variables. The spallation 
contamination remaining after the spallation cut is diffi- 
cult to estimate due to the large statistical uncertainties, 
but no evidence for remaining background was found. 
For the SK-I/III combined sample, we could see any ex- 
cess with a resolution of about 4 events; and for SK-II, 
with a resolution of about 2 events. 

SK-1 and SK-1II use the same spallation cut. Due to 
the differences in cathode coverage, SK-II required sepa- 
rate likelihoods and tuning. The lower energy threshold 
of the cut is the energy below which keeping the sam- 
ple free of any definite spallation contamination becomes 
too inefficient to be viable. Fig. [S] shows the half-life 
and maximum energy of spallation products in SK; as 
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FIG. 3. SK-I/III data with likelihood functions overlaid for 
single through-going muons. Top shows transverse distance, 
bottom shows longitudinal distance. 
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FIG. 4. Example of a dE/dx plot. The red line indicates 
where along the muon track the candidate was reconstructed. 
This example has particularly good correlation. 



spallation products decrease in energy, they tend to have 
longer half-lives, and are more difficult to tag. Taking 
into account energy resolution, events at 14 MeV such 



5 



o 1 ° 


„ 4 

Be 


ii 


Half-li 


: 16 n 


: 


1 


Ifi C 8 Li 
. 8 B • ' 


- 


-1 

10 


Li 

f S He 


- 


-2 

10 


12 B 

l2 Be' 

"B l2 N 
- ' ."() 


14 B 

"Li . 



6 8 10 12 14 16 18 20 22 

Energy (MeV) 



FIG. 5. Half-lives and event end-point energies of spalla- 
tion products expected to occur in pure water. Lower energy 
events trend to longer decay times. 



as 9 Li can reconstruct up to 16 MeV. The final energy 
threshold of the analysis, for which no evidence of spal- 
lation remains in the final sample, is 16 MeV positron 
energy (or 17.3 MeV D e energy) for SK-I/III, and 17.5 
MeV for SK-II. 

The cut need not be applied at all energies. The cut's 
upper energy bound is determined by the highest en- 
ergy spallation products, n Li and 14 B, each with an end 
point energy of 20.6 MeV and a half-life of around 0.01 
seconds. The SK-I/III spallation cut is applied up to 24 
MeV positron energy, and the SK-II cut up to 26 MeV, 
which is about 1.5 a of energy resolution away from 20.6 
MeV. 

The cut is split into two energy regions: a more strin- 
gent selection for the lower energies, which has many 
more spallation events and is most critical for the anal- 
ysis; and a looser selection at the higher energies, which 
has fewer spallation events, many of which are 11 Li and 
are easily removed considering its short half-life. Thus, 
for SK-I/III, the cut from 16-18 MeV is different than 
from 18-24 MeV, and in SK-II, the cut from 17.5-20 MeV 
is different than from 20-26 MeV. 

For all muons, a timing cut rejecting events with St < 
4 seconds (8 s) was implemented on the most energetic of 
muons, those over 400,000 (800,000) p.e., as these events 
are rare and almost always create spallation. 

For multiple muon type events, the secondary track is 
usually well fit, so the dE/dx profile of secondary tracks 
was calculated, and a separate secondary track likelihood 
tuned. Additional tracks do not have trustworthy dE/dx 
profiles, and a likelihood for these cases was constructed 
based only on At, £tranSj and the total amount of light 
deposited by the muon in the detector. 

Finally, for single through-going muons, ~ 1.5% of the 
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FIG. 6. Comparison of likelihood distributions of data sample 
and random sample for single through-going muons in SK- 
I/III. The data sample contains spallation events, while the 
random sample does not. The high likelihood values of the 
spallation excess show the effectiveness of the cut. The SK- 
I/III 16-18 MeV cut value of 3 is shown; everything to the 
right of the line is rejected. 



TABLE I. Spallation Cut Signal Efficiencies 



EFFICIENCY 


SK-I/III 


SK-II 


Low energy 
High energy 


81.8 % (16-18 MeV) 
90.8 % (18-24 MeV) 


76.2 %(17.5-20 MeV) 
88.2 %(20-26 MeV) 



time the muon fitter was unable to correctly fit the track. 
A dedicated alternate fitter was developed for these misfit 
events, which correctly fit approximately 75% of the mis- 
fit single through-going muons. For the remaining 25% 
of misfits, a 2 second detector veto was implemented to 
prevent spallation contamination. 

Fig. [5] demonstrates the power of this method to tag 
spallation events. Table fl] shows the cut signal efficiency. 
In the published 2003 SK analysis, in order to be spal- 
lation free the spallation cut was 37% inefficient from 18 
to 34 MeV. Now, for the same SK-I data, the cut is only 
9% inefficient from 18 to 24 MeV, and furthermore the 
16-18 MeV region is now opened up for study with an 
inefficiency of 18%, all spallation free. 



Solar angle cut 

Solar neutrinos elastically scatter off electrons, creat- 
ing recoil electrons that look like SN relic positrons. To 
discriminate such events, we use the angle between the 
reconstructed direction vector of the event and the direc- 
tion vector from where the sun was in the sky at the time 
of the event (0 Bun ). 

Although the direction of the recoil electron closely fol- 
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FIG. 7. SK-I/III solar MC integrated cos(6» 3U „) distributions 
by multiple Coulomb scattering goodness bin. 

lows that of the incoming solar neutrino (usually within 
~ 10 degrees), multiple Coulomb scattering (MCS) of the 
recoil electron smears out the peak in the cos(6> BU11 ) distri- 
bution at cos(6> BU11 ) = 1. This long tail in the distribution 
makes cutting solar neutrino events using cos(0 Hun ) ineffi- 
cient, especially since the end points for 8 B and hep solar 
neutrinos are right at the critical lower energy edge of the 
analysis, where inefficiency impacts the relic sensitivity 
the most. In order to maximize sensitivity, we used an 
energy dependent cut coupled with an estimate of MCS. 

MCS was estimated using a PMT hit by PMT hit 
Hough transform that reconstructs the direction and re- 
turns a goodness value corresponding to the amount of 
MCS. The correlation of this MCS goodness to the shape 
of the cos(# sun ) distribution can be seen in Fig. The 
candidate sample is split into various MCS goodness bins, 
and the cut value for each MCS goodness bin is separately 
tuned, in one MeV energy bins, by use of the following 
significance function: 

Significance = t/^/nS + ea (1) 

where e is the cut efficiency, S is the number of solar 
neutrino events, k is the reduction effectiveness of the cut 
(such that k x S is the number of solar neutrino events 
remaining after the cut is applied), and a represents the 
non-solar neutrino background. 

The solar neutrino spectrum is modeled using Monte 
Carlo (MC). The solar neutrino flux was normalized us- 
ing real SK data below energy threshold (14-16 MeV), 
with all cuts applied (except the one on cos(0 Hun )). To 
normalize, it was assumed that all events with cos(6* 5un ) 

< were not solar neutrino events, and that all non-solar 
remaining events had no cos(# aun ) dependence. Thus, 
simply subtracting the number of events with cos(# aun ) 

< from the number of events with cos(6* aun ) > allows 



TABLE II. SK-I/III solar cos(0 aun ) cut value (inefficiency) 



MCS Goodness 


16-17 MeV 


17-18 MeV 


18-19 MeV 


g < 0.4 
0.4 < g < 0.5 
0.5 < g < 0.6 

0.6 < g 


0.05 (7.2%) 
0.39 (9.7%) 
0.59 (6.7%) 
0.73 (2.4%) 


0.35 (5.0%) 
0.61 (6.4%) 
0.73 (4.5%) 
0.79 (2.0%) 


0.45 (4.3%) 
0.77 (3.8%) 
0.81 (3.2%) 
0.91 (1.0%) 


Total 


26.2% 


17.9% 


12.2% 



us to estimate the number of solar neutrino events at a 
particular energy, and then extrapolate the number of 
solar neutrino events at all energies. 

The dominant non-solar neutrino background was as- 
sumed to be decay electrons from atmospheric CC 
events, which is our largest background (see section III). 
For the purpose of tuning the solar cut, the decay electron 
spectrum was modeled using a sample of data that was 
identified as decay electrons by timing and spatial cor- 
relation to a preceding muon that had no OD trigger (a 
fully contained muon from a muon neutrino interaction). 
The normalization of the decay electron background was 
determined by fitting to the data. 

As with many of the cuts, SK-I and SK-III used one 
cut, and a separate version of the cut was tuned for SK-II. 
The SK-II cut criteria is simply the SK-I/III cut with the 
energy bins shifted by +6% to reflect the poorer energy 
resolution. The cut is MCS goodness bin dependent for E 
< 19 MeV (SK-I/III, < 20.14 MeV SK-II), while for 19- 
20 MeV (20.14-21.2 MeV SK-II) there is a simple, MCS 
goodness bin-independent cut that eliminates events with 
cos(# aun ) > 0.93. More cut details can be found in Table 
HIl Note that the cut values in the table mean that events 
with cos(# Bun ) > cut value are rejected. Also note that 
the inefficiencies are already weighted by the fraction of 
events in that MCS goodness bin, so that they can be 
simply added for the total. 

After the cut, no statistically significant evidence for 
remaining solar neutrino events exists, and the estimated 
number of remaining solar neutrino events in the com- 
bined SK-I/II/III final data sample is < 2. 



Incoming event cut 

The FV cut eliminates most of the radioactivity events 
that enter from detector edges, but not all. Instead of 
reducing the FV, it is more efficient to separate the re- 
maining radioactive background by position and direc- 
tion considerations. 

We use an effective distance parameter, d off which is 
defined identically as in the SK solar neutrino analysis 
(where this is referred to as the 'Gamma Ray Cut' |27j|). 
The event direction and position are reconstructed; then, 
starting at the event vertex, the distance projected back- 
wards along the event direction is calculated until a wall 
is hit. This distance is d eH . 
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FIG. 8. SK-I data distribution showing energy dependent 
incoming event cut. Events below the line are rejected. 



The radioactive backgrounds targeted by the incoming 
event cut are energy dependent, and much more prevalent 
at lower energies. However, some background exists at 
all energies. Thus, the cut has two pieces. A 300 cm d eff 
cut is applied at all energies, for all three SK phases. The 
lower energy part of the cut is tuned separately for each 
SK phase. 

For SK-I, sufficient statistics existed to allow an en- 
ergy dependent tuning of the cut, to maximize efficiency 
(see Fig. [S]). The shielding added to the PMTs at the 
start of SK-II to prevent another serious accident pro- 
duce radioactive backgrounds, and thus the SK-I tuning 
could not be used for other phases. Furthermore, the re- 
duced cathode coverage of SK-II makes it different. Each 
phase is tuned separately, but for SK-II and SK-III, the 
statistics were not sufficient for energy dependent tuning. 
Instead, a cut was implemented below 22 MeV, at d eft < 
500 cm for SK-II, and d eS < 450 cm for SK-III. 



triggered event within the 50 /is causes the candidate to 
be rejected. Since SLE events are much more frequent, a 
further vertex correlation requirement of 5 meters is re- 
quired in order to improve efficiency. When searching for 
post-activity, the 5 meter vertex correlation is required 
of all events, whether HE, LE, or SLE. 



Pion cut 

Some of the higher energy events in our sample (usu- 
ally > 30 MeV, reconstructed assuming event is an elec- 
tron) are charged pions from atmospheric neutrino inter- 
actions. One way to distinguish between charged pions 
and electrons is by looking at the 'fuzziness' of the re- 
sulting Cherenkov ring. This is because pions are quickly 
captured by an oxygen nucleus, while electrons tend to 
travel farther and MCS more. This MCS causes the 
electron direction to change as it travels, causing the 
Cherenkov ring to be more poorly defined. Pions are 
also heavier, and deflect less when scattering. 

Quantitatively, this difference can be used to construct 
a tool for discriminating pions and electrons. When the 
Cherenkov angle for an event is calculated, a distribution 
of three-hit combinations (or triplets) is formed. Each 
combination of three hit PMTs form a cone, and an open- 
ing angle. The distribution of these opening angles has a 
peak that is narrower for pions than for electrons. Using 
this logic, a pion ratio variable was constructed, where 
pion likelihood is defined as follows (using the triplet dis- 
tribution): 



Rpion — 



^triplets ± 3°from peak 
(#triplets ± 10° - #triplets ± 3 C 



(2) 



By looking at MC compared to data, the cut criterion 
was determined to be 0.58 for all SK run periods. 



Decay electron cut 

The SK electronics saves PMT signals in a 1.3 fxs win- 
dow. Sometimes, a muon decays into an electron quickly 
enough that the light from both particles are captured in 
the same event. These events are eliminated by searching 
for a multi-peak structure in the timing information of 
the event. 

More commonly, the muon and resultant decay elec- 
tron are in separate events. For each relic candidate, all 
SK data is searched up to 50 /is before the candidate (in 
case the candidate is a decay electron), and up to 50 /is 
after the event (in case the candidate is actually a low 
charge muon from an atmospheric neutrino interaction, 
or one of two decay electrons). 

When searching for pre-activity, finding any LE or HE 



Cherenkov angle 

Charged particles traveling in pure water at a veloc- 
ity close to the speed of light in vacuum emit Cherenkov 
light in a cone shape with an opening angle of about 42 
degrees. Heavier particles, such as muons, pions, and 
protons, may have a velocity that is smaller if their ener- 
gies are relatively small, producing Cherenkov light with 
an opening angle that is less than 42 degrees. Fig. [S] 
shows the Cherenkov angle distribution of SN relic MC 
(all models are very similar). The Cherenkov angle used 
is the peak of the distribution of three-hit combination 
opening angles. The signal region is determined to be 
between 38 and 50 degrees. Events at other Cherenkov 
angle values are not discarded, but instead they are kept 
to help extrapolate backgrounds into the signal region, 
as discussed in the following section. 
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FIG. 9. Expected Cherenkov angle distribution of SRN 
events. This is expected to be model independent; shown 
is Ando et al.'s LMA Q relic MC. 
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FIG. 10. Signal efficiency of new analysis compared to 2003 
study. An unnormalized LMA spectrum is overplotted to 
show the importance of the improvements. 



TABLE III. Cut Signal Efficiency (Systematic Error) 



Cut 


SK-I 


SK-II 


SK-III 


Noise Reduction 

Spall + Solar 
Incoming event 

Pion 
Cherenkov angle 
Other cuts 


99% (1%) 
88% (1%) 
98% (0.5%) 
98% (0.2%) 
95% (0.4%) 
98% (2%) 


99% (1%) 
87% (1.4%) 
95% (0.3%) 
97% (0.5%) 
88% (3%) 
98% (2%) 


99% (1%) 
89% (1%) 
96% (0.3%) 
98% (0.5%) 
94% (0.3%) 
98% (2%) 


Total 


78.5% (2.5%) 


69.2% (4.0%) 


76.7% (2.5%) 



Other cuts 



Occasionally events traveling through the OD do not 
deposit enough light to cause an OD trigger, but still 
cause a few OD PMT hits. OD hit clusters are compared 
in time and space to ID hits and the reconstructed event 
vertex, and if sufficient correlation exists, the event is 
rejected. 

The multi-ring cut uses ring counting software devel- 
oped for atmospheric neutrino and proton decay analyses 
^29] . If the software finds that more than one Cherenkov 
ring exists in the event, the angle between the two rings 
is calculated, as sometimes the fuzziness of low energy 
electron Cherenkov rings can cause the event to be mis- 
takenly identified as having more than one ring. False 
2-ring results reconstruct as having a small relative an- 
gle. A 2-ring result separated by more then 60 degrees is 
judged to be a real multi-ring event, and such events are 
rejected. 



B. Efficiency and systematic error 

Most cut efficiencies were calculated using relic MC. 
Each relic model has slightly different efficiencies; the 
ones listed in Table IIIII arc from the LMA model. The 
systematic error of most cuts was calculated by compar- 
ing LINAC MC to LINAC data, or by using the 
values established in the solar analyses when the cuts are 
similar (noise reduction, incoming event cut). The OD 
correlated cut uses a random sample of cosmic ray muon 
data for OD hit information for determining efficiency. 
The efficiency of the spallation cut is a function of posi- 
tion (efficiency is greater near the edges of the detector, 
where the candidate can be correlated to fewer muons), 
and the vertex distribution of relic MC is used to esti- 
mate the average efficiency. The error on this estimation 
is likely dominated by the statistics of our relic MC and 
spallation random sample, and is rounded up to be con- 
servative. The efficiency of the solar angle cut is mostly 
geometric and exact; the only source of systematic error 
comes from the ratio of events falling in each multiple 
Coulomb scattering goodness bin, as determined from 
solar MC. Values of the systematic errors, along with the 
cut efficiencies, can be seen in Table [LTTl 



III. REMAINING BACKGROUNDS 

A. Background Sources 

Even after all cuts, some backgrounds remain that 
must be modeled. The remaining background can be 
categorized into four main groupings. Atmospheric 
and v e charged current (CC) backgrounds were consid- 
ered in the 2003 SK published SRN paper, and now ad- 
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ditional atmospheric neutrino backgrounds are also con- 
sidered. Although remaining backgrounds other than the 
four listed below were found to exist (for example, from 
multiple and neutral pion production), their contribu- 
tions were small and their spectrum shapes well described 
by linear combinations of our four modeled backgrounds. 
Thus, any contribution from these sources will be ab- 
sorbed into the four backgrounds below during the fit. 



1: Atmospheric f M CC events 



This is the largest remaining background in our sam- 
ple. Atmospheric v^'s and P^'s interact in the water of 
the detector and create a muon via a charged current 
reaction. This muon is very low energy, often below 
Cherenkov threshold, in which case its decay electron 
cannot be removed like other decay electrons (by correla- 
tion to the preceding muon), since the muon is invisible. 
This background's energy spectrum is the well known 
Michel spectrum, slightly modified by resolution effects, 
and is quite different from the SN relic spectrum. 



C/5 




70 80 90 
Energy (MeV) 



FIG. 11. Spectra of the four remaining backgrounds in the 
signal Cherenkov angle region with all reduction cuts applied. 
The Vfj, CC channel is from decay electron data; the other 
three are from MC. All are scaled to the SK-I LMA best fit 
result. 



2: i/ e CC events 



This background also originates from atmospheric cos- 
mic ray interactions. Atmospheric ^ e 's and z/ e 's are in- 
distinguishable from SN relic i/ e 's on an individual basis. 
Their spectrum is quite different, however. 

3: Atmospheric v neutral current (NC) elastic events 

NC elastic events have an energy spectrum that rises 
sharply at our lower energy bound, similar to SN relics. 
Most are removed by Cherenkov angle reconstruction, 
but some still leak into our final sample and must be 
modeled. With the lowering of the energy threshold from 
18 MeV to 16 MeV, this background has become much 
more relevant. 



4: fj,/n production from atmospheric v 

This last category is a grouping of two different things, 
both heavier particles. First, NC reactions produce 
charged pions > 200 MeV/C, some of which survive the 
pion cut. These remaining events must be modeled. In- 
cluded with them, since the spectrum and Cherenkov 
angle distribution are relatively similar, are surviving 
muons above Cherenkov threshold. 



B. Background Modeling 

All our backgrounds are modeled using the SK MC, 
based on GEANT 3 and NEUT 30 . Hadronic interac- 



tions are simulated using a combination of CALOR and 
custom SK code. For the Michel spectrum (used for the 
CC), we measure the spectrum of decay electrons from 
cosmic ray muons. The spectra of the other three re- 
maining backgrounds are determined from the MC, as is 
the Cherenkov angle distribution of all four backgrounds. 
The SN relic signal itself is modeled separately, with its 
own MC. Many different models are considered, and dif- 
ferent results are calculated for each model. 

Fig. [TT] shows the energy spectrum of the four re- 
maining backgrounds taken from the MC. Fig. [T2l shows 
the Cherenkov angle of the four backgrounds, with no 
Cherenkov angle cut applied (from MC); it can be seen 
that the MC has a shape quite similar to the data. Also, 
while the CC backgrounds have a Cherenkov angle dis- 
tribution similar to that of SN relics (38-50 degrees) , the 
NC elastic background mostly reconstructs at high angles 
(where the Cherenkov angle reconstruction algorithm, 
which assumes a single ring, misreconstructs events with 
multiple gammas and roughly isotropic distributions of 
light), and the /i/tt leakage mostly reconstructs at lower 
angles (as expected from heavier, low energy particles). 
Because of these distributions, we have split the data into 
three regions: the signal region (38-50 degrees); and the 
two background regions, or 'sidebands', consisting of the 
low (20-38 degrees) /j,/tt region, and high (78-90 degrees) 
NC elastic region. The sidebands are used to normalize 
the NC elastic and /i/tt backgrounds in the signal region, 
where the SN relic signal would occur. 
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FIG. 12. Cherenkov angle of SK-I combined final data (all 
cuts except Cherenkov angle cut applied) overlaid with dis- 
tributions of the four remaining backgrounds from SK-I MC 
(same cuts applied). The division of Cherenkov angle regions 
is also pictured. 
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FIG. 13. Example likelihood curves from the LMA model. 
Vertical lines represent best fit and 90 % C.L. results for the 
combined fit. 



relationship, where r represents the number of SN relic 
events: 



IV. LIMIT EXTRACTION 



limit 90 poo 

£(r)dr/ / £(r)d> = 0.9 (4) 
ii Jo 



A. Likelihood fit 



The relic best fit and upper flux limit are determined 
by performing an unbinned maximum likelihood fit. A 
simultaneous fit is done in all three Cherenkov angle re- 
gions. For each of the Cherenkov angle regions the spec- 
trum of each of the five parameters (SN relic + 4 back- 
grounds) is parameterized into an analytical function, 
which is used as the PDF for that parameter. Each pos- 
sible reasonable combination of parameters is checked, 
and the best fit is the combination that maximizes the 
likelihood, which is determined by the function: 

Nevente 5 5 

£ II e-^U* x^CipdfjiEi) (3) 
»=i i=i 

where c is the coefficient for parameter j represent- 
ing the number of events for that channel present in all 
three Cherenkov angle regions combined, which modi- 
fies the PDF for that channel, pdf(E). The likelihood in 
each Cherenkov angle region is separately calculated, but 
maximized in conjunction. The likelihood is maximized 
for each SK phase separately. 

The 90% C.L. flux limit is extracted from the likelihood 
curve (maximum likelihood as a function of number of 
relic events). The likelihood curves of SK-I, II, and III 
are first multiplied together (see Fig. [TBI , then the 90% 
C.L. point limitgo is determined by the following simple 



B. Systematic error 

Four systematic errors are considered: 



1. Energy Scale and Resolution Uncertainty 

The energy scale and energy resolution systematic un- 
certainties are considered to be uncorrelated. They are 
separately calculated, then added in quadrature. The 
amount of the energy scale and energy resolution system- 
atic errors are based on the SK solar analyses [23, HH , 
then increased due to the broader energy range and 
higher energies, and to cover slight uncertainties in the 
MC. The energy resolution systematic uncertainty used 
in our study is 3%, while the energy scale uncertainties 
used are 2% for SK-I and SK-III, and 3% for SK-II. In- 
clusion of this systematic error was not found to have an 
appreciable effect on the results. 

To take the energy scale systematic uncertainty into 
account, the PDFs for four of the five channels (the 
CC was taken from decay electron data, and is thus not 
MC, and not susceptible to this error) were shifted by 
the amount of the uncertainty. To take into account the 
energy resolution systematic uncertainty, the energy res- 
olution function from reference [28| was applied to the 
likelihood functions to distort them. 
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Finally, the full uncertainty was incorporated into the 
likelihood calculation algorithm as follows. Let e be a 
variable representing the amount of spectral distortion, 
such that an e of 1 represents a lcr deviation in the en- 
ergy scale and resolution; and let £(e) be the likelihood 
evaluated with the likelihoods distorted by an amount e. 
Then the new final likelihood C is: 

C' = ^= e^/ 2 C(e)de (5) 

V ^ J-oo 

such that the originally calculated likelihood C is now 
simply the e = case. 

2. Atmospheric v Background Systematic Errors 

Three of the four remaining backgrounds are modeled 
by SK MC (the exception being the CC background, 
which is modeled by decay electron data). The SK MC 
is well verified above 100 MeV; however, in the lower en- 
ergy region of 16-100 MeV that is relevant to our study, 
the MC accuracy is less well studied. Little data exists to 
compare to the MC, especially for the NC elastic mode. 
To incorporate these concerns into the analysis, two sys- 
tematic errors on the v e CC and NC elastic backgrounds 
were included. These systematics were designed to cover 
any remaining uncertainties in the MC. As for the other 
two background channels, the CC comes from real 
data and therefore is not vulnerable to MC inaccuracies, 
and the fi/n channel is by far the least relevant of the 
four channels, and therefore the safest to neglect. 

For the NC elastic channel, the greatest potential for 
error lies in the relative normalization across Cherenkov 
angle regions. The amount of NC elastic found in the 
78-90 degree region determines the amount of the fit in 
the signal region. If the ratio of the amount of NC elastic 
in the high angle region compared to the signal region as 
determined by MC is off by some factor, this will tend 
to have a larger influence on the result than a distortion 
of the spectrum of the NC elastic in the signal region by 
the same factor. 

For SK-I, the NC elastic MC has 7.4% of events occur 
in the signal region, while 87% of events occur in the high 
angle region. Considered is a 100% change in the normal- 
ization of the number of NC elastic events in the signal 
region, with all changes to be correspondingly balanced 
by a change in the normalization of the high angle region. 
Thus, a +lcr effect would be for 14.8% of events to occur 
in the signal region, and 79.6% of events to occur in the 
high angle region, while a -lcr effect would be for no NC 
elastic events to occur in the signal region, and 94.4% 
of the events to occur in the high angle region. As we 
can't have less than events in the signal region, there is 
a physical boundary at -lcr. Thus, we applied the error 
asymmetrically from -1 to +3 a. 

For the v e CC channel, the relative normalization 
across Cherenkov angle regions is not a concern, as more 



than 99% of the events fall into the signal region, so in- 
stead a distortion of the spectrum in the signal region 
was considered. The v e CC PDF in the signal region was 
distorted as follows: 



PDF neu , = PDF oM x N(l + \l MeV ' ) (6) 

where E is the energy (MeV), e is the magnitude of 
the spectral distortion, and N is a constant such that the 
new PDF is normalized the same as the old. 

A -2cr effect would reduce the PDF to at 90 MeV. 
This was considered to be unphysical, and thus a -lcr 
effect was considered the physical lower bound, forcing 
us to apply this error in the same asymmetric manner as 
the NC elastic, from -1 to +3cr. 

We applied the background channel systematics using 
an integral method with a weighting function (Bayesian 
prior). However, as our error is applied over an asym- 
metric region (-1 to +3 a), we could not use a symmet- 
ric Gaussian (such as we use in the next section) with- 
out introducing a bias in our result. Instead, we used a 
weighted Gaussian that as closely as possible maintained 
the properties of the symmetric Gaussian. 

Specifically, in the case of a symmetric error (i.e., - 
oo to oo), a normal Gaussian would have the following 
properties: (1) expectation value (first moment) = 0, (2) 
variance (second moment) = a 2 . Our weighted Gaussian 
was constructed such that its properties (integrated over 
cr from -1 to 3) were the same as the symmetric case. 
This weighted Gaussian was normalized to 1 and used in 
both the v e CC and NC elastic cases, as their integration 
ranges were identical. 

The integral method of applying the systematic error 
is as follows. Let G(cr) be our weighted Gaussian func- 
tion described above. Let £(r, a) be the likelihood as 
a function of SN relic events and whichever systematic 
error is under consideration. Then, the likelihood after 
the application of the systematic error of either of the 
background channels is: 



£(r)new = J C(r,a)oi d G(a)da (7) 

Inclusion of the atmospheric v background systematics 
changed the result (towards a less stringent limit) by less 
than 6%. 



3. Energy Independent Efficiency Systematic Error 

The energy independent portion of the efficiency 
sytematic error is assumed uncorrelated. The system- 
atic error of each cut, as shown in Table IIII1 is added 
in quadrature, to get the total efficiency systematic error 
of the reduction. The total efficiency systematic error 
also includes the uncertainty in the fiducial volume cut, 
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TABLE IV. Total Efficiency Systematic Error 



Error Source 


SK-I 


SK-II 


SK-III 


Cut reduction 


3.1% 


4.4% 


3.1% 


Fiducial Volume 


1.3% 


1.1% 


1.0% 


Cross Section 


1.0% 


1.0% 


1.0% 


Livetime 


0.1% 


0.1% 


0.1% 


Total 


3.51% 


4.65% 


3.41% 



TABLE V. 90 % C.L. flux limit (y cm -2 s" 1 ), E v > 17.3 MeV 



Model 


SK-I 


SK-II 


SK-III 


All 


Predicted 


Gas Infall (97) 


<2.1 


<7.5 


<7.8 


<2.8 


0.3 


Chemical (97) 


<2.2 


<7.2 


<7.8 


<2.8 


0.6 


Heavy Metal (00) 


<2.2 


<7.4 


<7.8 


<2.8 


< 1.8 


LMA (03) 


<2.5 


<7.7 


<8.0 


<2.9 


1.7 


Failed SN (09) 


<2.4 


<8.0 


<8.4 


<3.0 


0.7 


6 MeV (09) 


<2.7 


<7.4 


<8.7 


<3.1 


1.5 



inverse beta cross section, and livetime calculation, as 
shown in Table [TV] 

The fiducial volume and livetime calculation system- 
atic errors are the same as, and taken from, the SK solar 
analyses [USl]. 

To apply this systematic error, the likelihood curve 
(likelihood as a function of number of SN relic events) 
is modified. The modification for each value of the like- 
lihood is as follows, where e is the efficiency, a is the 
efficiency systematic error, P(e) is a probability function 
in the shape of a Gaussian centered on eo with width eo<7, 
r is the number of relic events seen in the data, and R 
is the number of relic events actually interacting in the 
detector, such that r = eR: 



C'(R) = ( C(eR)P(e)ede 

Je=0 



(8) 



This error is applied last, after all other calculations are 
complete. Inclusion of this systematic error changed the 
result (making the limit less stringent) by a few percent. 



V. THE RESULT 

Figures [T4] [15] and [16] show the final sample data, 
overlaid with the best fit results for the LMA model. For 
SK-I, the relic best fit is negative and unphysical (for all 
models), so zero relic contribution is shown. However, 
for all models SK-II and SK-III find a positive indication 
(not significant) of a relic component. The 90% C.L. up- 
per flux limits differ for each model, and are summarized 
in Table IV] 




25 50 75 25 50 75 25 50 75 
20-38 degrees (MeV) 38-50 degrees (MeV) 78-90 degrees (MeV) 

FIG. 14. SK-I LMA best fit result. The relic best fit is nega- 
tive, so a relic fit of is shown. 




25 50 75 25 50 75 25 50 75 
20-38 degrees (MeV) 38-50 degrees (MeV) 78-90 degrees (MeV) 

FIG. 15. SK-II LMA best fit result. The relic best fit is 3.05 
events per year interacting in the detector (before reduction 
efficiencies). 




25 50 75 25 50 75 25 50 75 
20-38 degrees (MeV) 38-50 degrees (MeV) 78-90 degrees (MeV) 

FIG. 16. SK-III LMA best fit result. The relic best fit is 6.9 
events per year interacting in the detector (before reduction 
efficiencies). 



A. Comparison to previous study 

The new analysis has a lower energy threshold than the 
2003 study, so the results are not immediately compara- 
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FIG. 17. True positron spectra in SK for each neutrino tem- 
perature, from 3 MeV to 8 MeV in 0.5 MeV steps (SN P e 
luminosity of 5xl0 52 ergs assumed). 
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ble. Furthermore, the previous study presented its result 
in a model independent fashion due to the similarity of 
the results. Although we verified that with an 18 MeV 
energy threshold the SK-I results were similar enough 
to justify a model-independent approach, the lower en- 
ergy threshold and new data required model dependent 
results. For comparison purposes, we consider the LMA 
model. The limit is 2.9 v cm" 2 s _1 > 16 MeV (positron 
energy), which is equivalent to 2.0 v cm~ 2 s _1 > 18 MeV 
positron energy. It is notable that this result is less strin- 
gent than the 2003 result of 1.2 v cm~ 2 s _1 positron 
energy > 18 MeV. There are multiple reasons for this. 

First, a th order approximation of the inverse beta 
cross section was then used. Now, the full cross section 
from 25] is used. This raises the limit by about 8%. 
If events with post-activity are also removed, the old- 
style analysis limit becomes 1.35 cm s . Furthermore, 
the binned \ 2 method used assumed Gaussian statistics, 
while Poissonian statistics are more appropriate consider- 
ing the low statistics. This alone would change the limit 
from 1.2 to 1.7 cm -2 s — 1 . When all these corrections are 
combined, the original analysis result of 1.2 v cm -2 s _1 
instead becomes 1.9 v cm -2 s _1 . 

With our improved analysis, if we neglect atmospheric 
v background systematics (which were not fully included 
in the 2003 study), the SK-I only LMA result is 1.6 v 
cm -2 s _1 (> 18 MeV positron energy), which is more 
stringent than the published analysis with these correc- 
tions. However, the SK-II and SK-III data shows a hint 
of a signal, which causes the limit to become less strin- 
gent when all the data are combined, for the final LMA 
result (with all systematics) of 2.0 v cm~ 2 s _1 > 18 MeV 
positron energy, or 2.9 v cm~ 2 
energy. 



> 16 MeV positron 



B. Typical SN v emission limit 

Most of the elements involved in a comprehensive pre- 
diction of the SRN flux are now fairly well known (32| 



FIG. 18. Results plotted as an exclusion contour in SN neu- 
trino luminosity vs. neutrino temperature parameter space. 
The 1MB and Kamiokande allowed areas for 1987A data are 
shown (originally from [3j|) along with our new 90% C.L. re- 
sult. The dashed line shows the individual 90% C.L. results 
of each temperature considered separately, which is not a true 
2-D exclusion contour. Results are in the form of Fig. 6 from 

m. 



(e.g. initial mass functions, cosmic star formation his- 
tory, Hubble expansion, etc.), and thus we can parame- 
terize typical supernova neutrino emission using two ef- 
fective parameters [l3l [33| : D e luminosity from a typical 
supernova, and average emitted v e energy. 

To calculate this result, SRN MC samples were created 
with a Fermi-Dirac (FD) emission spectrum at multiple 
temperatures (see Fig. [T7]). The Ando et al. model @ 
corresponds to a FD spectrum of close to 6 MeV and a 
u e luminosity of 7xl0 45 J. MC was created from 2.5 to 8 
MeV in 0.5 MeV steps. PDFs were made of these spectra, 
and the final data sample was fit, as previously described, 
with these new likelihoods in the place of the other SRN 
models. The following constants were assumed: J7 m = 
0.3, fl A = 0.7, and c/H = 4228 Mpc. The result is 
independent of the initial mass function used. 

Our result is shown in Fig. [T5] and Fig. Q1JJ Also 
plotted for convenience is the 1-D line representing the 
results of each neutrino temperature analysed separately 
(the numerical values are given in Table IVI|) . 



VI. CONCLUSION AND FUTURE 

In summary, we have made large improvements to the 
2003 SK relic analysis. The cuts have been optimized, 
greatly increasing efficiency, up to around ~ 75%, and 
lowering the energy threshold down to 16 MeV. New re- 
maining backgrounds are now considered and modeled, 
and we have almost twice as much data in our final 
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3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 

T in MeV 



FIG. 19. Exclusion contour plotted in a parameter space of 
SRN event rate vs. neutrino temperature. The red contour 
shows our 90% C.L. result. The dashed line shows the indi- 
vidual 90% C.L. results of each temperature considered sep- 
arately, which is not a true 2-D exclusion contour. CGI is 
Cosmic Gas Infall model, HMA is Heavy Metal Abundance 
model, CE is Chemical Evolution model, LMA is Large Mix- 
ing Angle model, FS is Failed Supernova model, and the 6 and 
4 MeV cases are from [131 ] . For the 4 and 6 MeV cases a total 
uncertainty is provided and shown, and the HMA model gives 
a range which is shown. Other models have no given range or 
uncertainty and are represented by a star. 



TABLE VI. Summary of limit values. 



u e Temperature (MeV) 


Energy (xl0 5-i ergs) 


3.0 


9.0 


3.5 


4.8 


4.0 


2.9 


4.5 


2.1 


5.0 


1.5 


5.5 


1.3 


6.0 


1.1 


6.5 


0.92 


7.0 


0.82 


7.5 


0.75 


8.0 


0.70 



sample. A new maximum likelihood fit with multiple 
Cherenkov angle regions is utilized to extract a set of 
model-dependent upper flux limits. Further, we have ex- 
pressed our results in a model-independent fashion using 
two effective parameters to model typical SN neutrino 
emission. 

Although our limit is tantalizingly close to the best the- 
oretical predictions, no signal has so far been detected. 
While SK is expected to continue data taking for many 
more years, future sensitivity improvements will be slow, 
since SK-I/II/IIFs exposure is already 176 kt-years, and 
the analysis is now highly optimized. 50 kt-years of SK- 
IV data already exists, and the new SK-IV electronics 
structure may allow for some further background reduc- 
tion, as decay electrons from atmospheric CC interac- 
tions can now be tagged by detection of prompt gamma 
emissions with higher efficiency. Further improvement 
would be possible with new methods, such as the doping 
of SK water with gadolinium, which could lower the en- 
ergy threshold and backgrounds dramatically and could 
allow detection of the SRN signal within five years [34| ■ 
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